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We consider evolutions of linear fluctuations as the background Friedmann world model goes from 
contracting to expanding phases through smooth and non-singular bouncing phases. As long as 
the gravity dominates over the pressure gradient in the perturbation equation the growing-mode 
in the expanding phase is characterized by a conserved amplitude, we call it a C-mode. In the 
spherical geometry with a pressureless medium, we show that there exists a special gauge-invariant 
combination $ which stays constant throughout the evolution from the big-bang to the big-crunch 
with the same value even after the bounce: it characterizes the coefficient of the C-mode. We 
show this result by using a bounce model where the pressure gradient term is negligible during the 
bounce; this requires additional presence of an exotic matter. In such a bounce, even in more general 
situations of the equation of states before and after the bounce, the C-mode in the expanding phase 
is affected only by the C-mode in the contracting phase, thus the growing mode in the contracting 
phase decays away as the world model enters expanding phase. In the case the background curvature 
has significant role during the bounce, the pressure gradient term becomes important and we cannot 
trace C-mode in the expanding phase to the one before the bounce. In such situations, perturbations 
in a fluid bounce model show exponential instability, whereas the ones in a scalar field bounce model 
show oscillatory behaviors. 

PACS numbers: 04.20.Dw, 98.80-k, 98.80.Cq, 98.80.Hw 


I. INTRODUCTION 

The collapsing and bouncing phases of the FLRW 
(Friedmann-Lemaitre-Robertson-Walker) world models 
we are considering are the possible ones in our past, 
before the big-bang. The same physics, however, could 
work in the possible case in future as well. 

The re-expansion of a positive curvature Friedmann 
world model which is destined to collapse, or cyclic rep¬ 
etition of the process with such a bounce was proposed 
as early as in the 1930’s [[j] ||. Specific realizations of the 
bounce and the conditions required to have the FLRW 
world model from a bounce were studied in The 

singularity-free cosmologies are possible as we give up 
the strong energy condition which is often possible with 
quantum corrections [Q. Recently, the big-bang world 
model preceded by a collapsing phase attracted renewed 
attention in the context of the brane cosmology ||. 

In this paper we analyse the evolution of scalar-type 
curvature (often called the adiabatic) fluctuations as the 
background world model goes through a smooth and non¬ 
singular bounce which connects the contracting and the 
expanding phases. We will assume the classical General 
Relativity is valid as the correct gravity theory through¬ 
out the evolution, and also consider scales where the lin¬ 
ear approximation is valid. 

In §[n] we review the cosmological perturbation the¬ 
ory needed for our analyses in later sections. In §111 
we present the large-scale evolutions of various curvature 
perturbations near singularity. In § IV we analyse the 


evolution using exact solutions in a pressureless situation. 
In §[y| we show the evolution of perturbation through a 
bounce model using three different bounce models. § VI 
presents a summary, implications of our work, and dis¬ 
cussions on related works. We set c = 1. 


II. COSMOLOGICAL PERTURBATIONS 


Our metric convention is 

ds 2 = —a 2 ( 1 + 2ot)drf — 2a 2 (/3 a + B^)dr]dx c 


+a z 


9 a 0 (1 + + 27,0,1/3 + 2G ( 


(v) 

a\/3) 


+2 C, 


dx a dx/. 


( 1 ) 


The perturbed order variables a , (3, ip and 7 are 
scalar-type perturbations; the transverse B^ 1 ' 1 and Ca ’■* 
are vector-type perturbations (rotation); a transverse- 
tracefree is tensor-type perturbation (gravitational 
wave). The energy-momentum tensor is 

To = ~{p + 6p), 

T ° = (h + p) -(iA>, a + ^ v) , 

T$ = (p+5p)8$+rf, ( 2 ) 

where the tracefree is the anisotropic stress. 

The trace and tracefree parts of extrinsic curvature 
(equivalently, the expansion 6 and the shear o a b of the 
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normal frame vector field), and the intrinsic scalar cur¬ 
vature of the constant-time spacelike hypersurface 
are, see eqs. (C3,C14) of || and eqs. (A6,A7) of Q, 


The scalar-type perturbation of a fluid with vanishing 
anisotropic stress in the Einstein’s gravity is described 
by 


(9 = 3 H-k, 

= x, a \ft - \$Ax + a4, (I|/3) + a2 CcS> 

R {h) = ^ [ 6 K - 4(A + 3 K)<p \, (3) 


where^j 


X = a((3 + cry), n = 3(Ha 
*£> = BW + aC£\ 



( 4 ) 


We have H = -, K the normalized background three- 
space curvature, and an overdot indicates time deriva¬ 
tive based on t, dt = adr/. Thus, k, X- and p are the 
perturbed expansion, the scalar-type shear, and the per¬ 
turbed three-space scalar-curvature of the normal hyper¬ 
surface, respectively. and C^g give the vector- and 
tensor-type contributions to the shear tensor. 

C^g , and Va' 1 are gauge-invariant, a, ip, x, k, v, 
Sp and Sp are spatially gauge-invariant but depend on 
the temporal gauge condition, i.e, depend on the spatial 
hypersurface (time slicing) choice |J. Setting any one 
of these temporally gauge dependent variables equal to 
zero corresponds to a fundamental gauge condition; ex¬ 
cept for the synchronous gauge (a = 0 ) each of the other 
conditions fixes the temporal gauge degree of freedom 
completely, and any variable in such a gauge condition 
uniquely corresponds to a gauge-invariant combination 
(of the variable and the variable used in the gauge con¬ 
dition) (|§] . 

Equations describing the evolution of a spatially ho¬ 
mogeneous and isotropic FLRW world model are: 

r2 = 87tG^ _ K_ + A _3H(p + p). ( 5 ) 

3 a z 3 


The A term can be also considered as an ideal fluid with 
Pa = — Pa = A/( 87 tG); in such cases, we have p = p m + 
Pa and p = p m + pa- In the case of a minimally coupled 
scalar field (</>) we have p = \(/> 2 + V and p = — V, 

with an equation of motion <j>+ 3H(f> + V^ = 0. A can be 
included as Va = A/( 87 rG). 

We consider the general scalar-type perturbation. It is 
convenient to introduce 10 |l 2 


= <Pv - 


K/a 2 


47 iG(p + p) 


Px- 


( 6 ) 


* Indices of d a g (E a g and -Ha/ 3 , later) are based on the 
spacetime metric g a b- All the other Greek indices are based 
on the g^. 


* H2 ( a -\ 


( 7 ) 

47 xG(p+p)a \H^ k ) 
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• He 2 k 2 

H 

( 8 ) 

4:TrG(p + p) a 2 

1 

M+P 

k 2 - 3 K 

2 Px—inGpdy, 


(9) 


where w = p/p and c 2 = p/p. p v = <p — ( aH/k)v , 
<p x = tp — H\ and S v = 5 + 3(aH/k)(l + w)v are gauge- 
invariant combinations ip v is the same as <p in the 

comoving gauge (v = 0 ), and p x is the same as ip in the 
zero-shear gauge (x = 0 ), etc. 

We emphasize that results in this section are valid con¬ 
sidering general K, A, and time-varying equation of state 
p(p). In the case of a minimally coupled scalar field, we 
have additional nonvanishing entropic perturbation (the 
isotropic stress), e = (1 — c 2 )5p v . Its effect can be covered 
by changingCg in eq. (@) to c\ = 1 — 3(1 — c 2 )K/k 2 , [ fl3|| ; 
as eqs. (@-@) are a complete set for single component, 
this prescription applies always in the single scalar field 
case. It is convenient to have 


W5y+ 3(ITW ,| ’ + g l2 <10) 

We can show in general ||l3[ 


Pn = Ps 


/ 


k 2 - 3 K 
127 tG(p + p)a 2 


( 11 ) 


where ip K = ip + Hn/(3H — k 2 /a 2 ). In the notation of 
Bardeen || we have 

dv — Cm, Px = 4*//, Pv — 0m, Pk = 4>h- (12) 

We have ips = C, in [@, and ip v = 1Z in Jh|; ip v was also 
originally introduced by Lukash as in Jl5| . From 
eqs. (@@) we can derive equations in closed forms: 

$ + (c 2 fc 2 /a 2 — x/x) $ = 0, (13) 

P x + (c 2 s k 2 / a 2 - y/y) <p x = 0, (14) 

where 

$ = <p x = (ay/H)p x , 

V = H /V(P + P) a = ( a/c 8 )x~ l . (15) 


Equations using the conformal time were presented in 

©■ 

In the large-scale limit (meaning c 2 k 2 /a 2 term neg¬ 
ligible, thus gravity dominates over pressure) we have 
general solutions |ll] Q : 
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t 


$(fc, t) = C{k) — d(k) 


4t tG 

t 


dt/x 2 , 


(16) 


ip x (k,t) = 4:TrGC(k)— f dt/y 2 -\-d{k) —, (17) 

CL J CL 


where C and d are the two spatially dependent integra¬ 
tion constants: we call these the C-mode and the d-mode, 
respectively. Solutions for 6 V , tps and tp K follow from eqs. 
(pHll|). Notice that the d-mode of $ is higher-order in the 
large-scale expansion compared with the d-mode of ip x . 
In a pressureless medium, the above solutions are exact, 
and we have $ = C | [10|. In fact, for such a medium, 
instead of eq. (13), eq. (||) gives $ = 0. 

In order to use the large-scale solutions in eqs. ( [hsUi?] ) 
it is important to check whether we could ignore the 
c 2 k 2 /a 2 term during the evolution. The large-scale con¬ 
dition implies (pressure/gravity) -C 1 where 


pressure c 2 k 2 /a 2 c 2 k 2 /a 2 

gravity x/x ’ y/y 


(18) 


In a positive curvature (spherical) model the wave num¬ 
ber varies as k = W (n 2 — \)K where n = 1,2,3,...; n = 
1,2 are known to be unphysical [I?li9j. In a negative cur¬ 
vature (hyperbolic) model, k > \J\K\ and 0 < k < \J\K\ 
correspond to the subcurvature and the supercurvature 
scales, respectively |l8). In a zero curvature (flat) model 
we have k > 0. 

The following two variables are continuous under a 
sudden jump of equation of state jlij: 

ip x (or <5„), (p s . (19) 


These joining variables work for general K , A, and p(y) 
in the general scale. This applies for the perfect fluids, 
and for the cases involving scalar fields, see (T^] . For 
the background, a and a should be continuous at the 
transition. Consider two phases / and II with different 
equation of states, wi and wn , making a transition at t\. 
Assuming a flat background, in the large-scale limit, by 
matching <p x and ips in eqs. (|l^,|l^.^o|) we can see that 
to the leading order in the large-scale expansion we have 


C n = C u 

du = dj + 47 tGOi 
f* 1 a(/i + p) 


fl a{n+p ) 
H 2 


dt 


H 2 


dt 


ii. 


( 20 ) 


Thus, to the leading order in the large-scale expansion 
the C-mode of $ remains the same, whereas the d-mode 
of is affected by the transition and also the previous 
history of the d- and C-modes |[u|. Applications were 
made in Q. 

Ignoring the anisotropic stress (7rS) and assuming K = 
0, equations for the rotation and the gravitational wave 
become (|] 


[a i {p + p)v^ ) ] =0; = 16 tt G(p+p)v£’ ) , (21) 

+ ( k 2 /a 2 - z/z) C® = 0, (22) 

where C^ = zC^p and z = a 3 / 2 . The equation for C^ 
satisfies the same equation as $ in eq. (0) with x oc a 3 / 2 . 
Thus, in general scale for and in the large-scale limit 
for we have general solutions 

*W(M) = d Q (fc)4, ( 23 ) 

a z 

Cg(/c,t) = c a0 (k) - d a0 (k) J dt/a 3 . (24) 

Similarly as the C-mode of $ the amplitude of c a( 3 -mode 
simply stays constant. 


III. LARGE-SCALE EVOLUTIONS OF 
CURVATURE PERTURBATIONS 

The general large-scale solutions for scalar- and 
tensor-type perturbations, and the general solutions 
for the vector-type perturbation are presented in eqs. 
dj - BHUg] ) . We call the solutions with C and c a0 the 
C-modes, and the solutions with d, d a and d a0 the d- 
modes. The vector-type perturbation has no C-mode. In 
the expanding phase C-modes are relatively growing solu¬ 
tions whereas d-modes decay, thus transient, in time. In a 
contracting phase, however, the opposite is the case with 
the d-modes often diverging as the background model 
approaches the singularity. 

In this section we assume near flat background. With 
a constant w we have c 2 = w, and a oc |f| 3 ( 1 +«") for 
— 1 < w < 1. In a medium with w > — as we ap¬ 
proach the singularity ^(oc 13q=35F) becomes negligible 
for any given scale with fc, thus the large-scale conditions 
are well satisfied. In such a case, during dynamical time- 
scale of the background evolution “light can travel only 
a small fraction of a wavelength” |9[, thus the scale be¬ 
comes super-horizon scale; Bardeen ||] called it (^ ~ 1) 
an “effective particle horizon”^]. In general, we consider 


f The effective particle horizon is the same as the “Hubble 
sphere” studied in [^l|, and closely resembles the “^-surface” 
introduced in [^2|. The global concepts like particle- and 
event-horizons are not suitable to describe the local dynami¬ 
cally reachable ranges. Studies in |^I|,^2| show that the Hub¬ 
ble sphere and the z-surface are more suitable to describe the 
concepts like ‘scales becomming super-horizon size during the 
inflation era’. Similarly, these are suitable to describe the 
same physics during contracting phase with w > — |. Under 
this situation we can show that to an observer in the contract- 
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w > 0: “a single-component treatment of the matter is 
inappropriate when the net pressure is negative” jl|. 

In this case, from eqs. (16 |l7| , |24| ) we have the C-modes 
remain constant in time 


‘Px = 


3 + 3w 
5 + 3w 


C, f v — C, 


— Ca 
° a0 ~ C «/3' 


(25) 


Thus, in an expanding medium the perturbation evo¬ 
lutions in the super-effective-particle-horizon are kine- 
maticr and are characterized by the conserved quantities, 
see [g4|. 

Meanwhile, the d-modes behave as 


Vx a i*r 


5 + 3u. 
‘ 3+311 




c *l «i*r 


i —ii 
■ l+ii 




( V ) 


\t\- 


3(l+») 


(26) 


for — 1 < w < 1. For w = 1 we havef^], from Table 4 of 

Vk, Va, VS, <Pv, OC In |f |, ip x oc |t| _4/3 . (28) 

Thus, even for w = 1, ip diverges in all gauge condi¬ 
tions considered. <p in the zero-shear gauge diverges more 
strongly compared with the ones in the other gauge con¬ 
ditions. The strong divergence in the zero-shear gauge is 
known to be due to the strong curvature of the hypersur¬ 
face (temporal gauge condition) f)]. ip is set to zero in 
the uniform-curvature gauge. For the C-modes we have 
eq. (|H) and tp K = tp a = tps = ip v = C. 


B. Extrinsic curvature 


for — 1 < w < 1. For w = 1, the above solution is valid 
for <£> x (oc |f| -4 / 3 ) and 4/^ (oc |t| 2 / 3 ), whereas we have 
ip v , oc ln|t| instead. 

In the case of constant w, a complete set of solutions of 
the scalar-type perturbation in six different fundamental 
gauge conditions is presented in Tables 2-5 of [|s|- Al¬ 
though the solutions in [H| are presented in the context 
of the expanding phase, by changing the time variables 
to their absolute values with the singularity at \t\ = 0, 
the same solutions apply in the contracting phase as well. 


A. Intrinsic curvature 

ip is a dimensionless measure of the (intrinsic) curva¬ 
ture perturbation of the hypersurface (temporal gauge 
condition) we choose. Thus, its value depends on the cho¬ 
sen hypersurface (temporal gauge condition). From eq. 
(A6) of § we notice that gives a dimensionless con¬ 
tribution to the tensor-type intrinsic curvature pertur¬ 
bation, and the vector-type perturbation does not con¬ 
tribute to the curvature perturbation, see also eq. (C14) 
in j|]. Table 2 of |2f| shows that for the d-modcj] 

Vn, Von VS, Vv, C y a ’p OC \t\ !+“, Vx a 1*1 3+3 ”> ( 27 ) 


Although we have no shear in the background model, 
the perturbed scalar-type shear is still a gauge dependent 
quantity. A dimensionless measure of the shear variable 
(the shear divided by the background expansion rate) 
becomes 


<7 

H 


f k 2 

\a 2 H 2 


H X , 


aH 


v&bO 

^ a j 


1 A(t) 


(29) 


for three perturbation types; a = yjaapa a P /2. The d- 
modes of all-types of perturbations show the same tem¬ 
poral behavior 


(30) 

for — 1 < w < 1, thus, with no logarithmic divergences 
for w = 1 case. The result for vector- and tensor-type 
perturbations follow from eqs. (|23|,HJ), and the one for 
scalar-type perturbation follows from the Tables 2 and 4 
of thus, the solutions apply to the gauge conditions 
considered (we set y = 0 in the zero-shear gauge). We 
can check that the C-modes contributions to a/H are all 
regular near the singularity for — 1 < w < 1. The be¬ 
havior of n/H, a dimensionless measure of the perturbed 
trace part of the extrinsic curvature, varies widely de¬ 
pending on the gauge conditions, see Tables 2 and 4 of 


ing phase, an object separated by a given comoving distance 
appears more blue-shifted as time goes on, i.e, the object is 
effectively receding from the observer. 

* This kinematic nature of the evolution is reflected in the 
conserved behaviors of linear perturbations in the expanding 
phase |^|^j. What is conserved is the the amplitude of the C- 
mode of curvature perturbation in diverse gauge conditions. 
In the expanding phase the situation is also well represented 
by the separated Friedmann picture of the perturbed Fried¬ 
mann world model pioneered (up to our knowledge) by Har¬ 
rison in Q. 

§ Although ip a , the i p in the synchronous gauge (a = 0), is 
not gauge-invariant, we are considering the physical solutions. 


C. Weyl curvature 

Durrer has informed us another useful measure of the 
spacetime fluctuation which behaves regularly for w = 1, 
the Weyl curvature C a bcd ■ The Wey 1-curvature (the con¬ 
formal tensor) vanishes in the FLRW background geom¬ 
etry, and is naturally gauge-invariant. The Weyl tensor 


** One other special case occurs for w = 0 where we have 
no d-modes for ip v and ip a , see Table 5 of [pH|. 
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can be covariantly decomposed to the electric E a s and 
magnetic H a b parts, [2^|. Using eq. (C9) in ||, see also 
eqs. (2.26,2.27) of we can show 


E 

R 


( k 2 


a 2 H 2 


09 — AW — <7^ A 

^ aH 2 ^ a ’ 


(31) 


Equations (p^|,pUp|-pd]) give exact solutions: 


< 1 > = C, 

(37) 

V>X - l P+ C + Vdd -,291 

k z — 3 1 — cos 77 

(38) 

PS = c + (k 2 / 9)(1 - cos 77 )^,, 

(39) 


where E = ^E a PE a p /2 and i? is the scalar-curvature 
(~ 77 2 ). Thus, the d-modes of all-types of perturbations 
behave exactly like a/H(oc |t| _T +"), thus behave regu¬ 
larly for w = 1. H a j, only contributes to the vector- and 
tensor-type perturbations, and we can also show that the 
d-modes behave as jH a ^H a p/2/R oc thus 

behave more regularly for w = 1. This regular behavior 
of the Weyl curvature at singularity for w = 1 case was 
used to argue the validity of perturbation theory in such 
a situation, see around eq. (5.20) of J|i|; however, see 
our discussion below eq. (|6§|). 

IV. EXACT EVOLUTION IN A PRESSURELESS 
CASE 

For I\ > 0, A = 0 and p = 0, eq. (|5|) gives a cycloid 

HI 

a = Cm (1 - cos 77 ), t = Cm (t? - sin 77 ), (32) 

where c m = (47rG/3)/ro 3 and dr] = dt/a. K is normalized 
to unity, thus 0 < 77 < 27r. 

For p = 0, eqs. (|],|b4[) give 

[a 2 H 2 (5 V /H)'Y /(a 2 H) 

= 5 V + 2 HS V + AnGpSy = 0, (33) 


and similarly for p K \ d = d/c ^ is dimensionless. For 
?7 <C 1 we have p x = |C and ps = p K = C for the 
C-modes. From eq. (||) we have 4) = 0 exactly for a 
pressureless fluid considering general K and A. $ has 
only the C-mode (it is identified as C), and no d-mode. 
Evolutions of p x and $ are shown in Fig. [j]. Notice that 
both p + and Pd diverge as the model approaches the big 
crunch singularity. 



v 


which coincides with the density perturbation equation in 
the synchronous gauge p~7| ], or in the Newtonian context 
p0| ; S v is the energy density perturbation in the comov¬ 
ing gauge, [[u]]. Assuming A = 0, the two independent 
exact solutions for p x oc S v /a in eq. © are Q: 


H 

1 

a 


dt 


sin 77 


a 2 H 2 

= ¥>+fa)/3, 


+ 


5 + cos 77 


H 


sm ?7 


m a (l-cos ??) 3 
In asymptotes we have 


(1 — cos r]) 3 (1 — cos r]) 2 

= <Pd(v)- (34) 


FIG. 1. Right: Evolutions of p+ (blue, line), p~ (cyan, 
dashed line), ps (blue, dotted line), and $ (black, horizontal 
line). 

Let us consider a scenario where the big crunch is suc¬ 
ceeded to an expanding phase: thus we have two phases 
V < Vi (phase I) and 77 > 771 (phase II) with 771 = 0. 
For p x we can take two out of three forms of solutions 
(p+, p~, Pd) as the general solutions in either phase. Al¬ 
though d could be discontinuous at the transition reach¬ 
ing the singularity let us try matching <p x and ps directly 
at 771 ; afterall we assume a nonsingular bounce near 771 , 
see later. Using eqs. (^,^9|,^5|) we can show 


<P+W — 3/5, ip d (ii) — 8 / 77 s , (77 <C 1), 

<P+(v)/ (487r) - (8/fj 5 ) ~ -<p d (i 7 ), (?7 < 1), (35) 

where fj = 2n — 77 . We have 

P+(2tt - 77 ) = ip+{rj) + 187079 ^( 77 ) = ( 77 ), (36) 

where p- shows time inverted evolution of <p+, |f33|. 


Cji = Cj , dn = dj — 187tCj. 


(40) 


This implies that the value of $ variable is conserved even 
through the bounce. In terms of the general solutions we 
notice the following. Using eq. (36) we can decompose 
p+ in eq. (^ 8 |) to p- and pd- Near bounce, although 
p[_ ~ | is negligible compared with p ! + ~ 187 t( 8 /|t 7 | 5 ) 
and p d ~ — 8 /I 77 I 5 , we can show that it is this constant 
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mode of <p{_ which feeds the C-mode after the bounce. 
Thus, in the collapsing phase it is appropriate to write 
eq. m as 

p x = ip-C + ipd{d — 187rC). (41) 

Therefore, if such a bounce is allowed, we have shown 
that < p{_ feeds the growing mode <p+ in the expanding 
phase. The apparent growing (diverging) solutions or 
p T d only feed the ip^ 1 or <pU which are the decaying mode 
in expanding phase. The time-scale of a cycle is encoded 
in c m of eq. @) which can affect only the decaying 
solution in the expanding phase. The value of <F, which 
is C, is not affected by the different duration of each 
cycle. 

Notice, however, that if we strictly consider the sin¬ 
gular and cuspy bounce at rji implied by eq. (^) we 
have a discontinuous, which forbids us from relying on 
the matching conditions. We have assumed such a singu¬ 
lar bounce can be regarded as a limiting case of a smooth 
and non-singular bounce; a concrete example will be con¬ 
sidered in the next section. We note that the curvature 
term has negligible role near the big crunch/bang. We 
also have assumed the linearity of the fluctuations in¬ 
volved. 


V. THROUGH THE BOUNCE 


Assume two expanding phases / and II with equation 
of states wi and wij. In near flat situation we have 


a(t) = o 0 i{t - U) 3<1+ ” 


(42) 


where i = I, II. The coefficients should be determined 
by matching a and a at the transitions t \. In phases / 
and II the large-scale solutions in eqs. (|1^,[L7 ]|lo|) give 
for — 1 < Wi < 1: 


$ = C + k' 


Wi 


1 


91-w? a 3 H 


d, 


3 + 3 wi , H 
<Px = c , o... C + — + 


5 + 3 Wi 


ips = C + k 2 - 


1 


1 


9 1 — Wi a 3 H 


(43) 

(44) 

(45) 


where C and d in the phase I should be regarded as Cj 
and di, and similarly for the phase II. For w t = 1, 
d-mode of $ (and part of tpg) contains In (t — ti) term in¬ 
stead of (1 —Wi) -1 . Equations (| 11) give <p v and ip K . As 
long as we have a and a continuous through a transition 
from phase / to phase II at t\ we can use our joining 
variables in eq. (|l^). Examples are the radiation-matter 
transition and the inflation-radiation transition. Using 
eqs. ( [i~i| . [f5|) , to the leading order in the large-scale ex¬ 
pansion we have 


C n = C7, 

dji = di + 


6 (wi - wjj ) ai 


C r . 


This is consistent with the result derived in eq. (17) of 
& Similar results hold for two contracting phases as 
well. 

In the case of transition from the contracting to ex¬ 
panding phases, however, a can be discontinuous at the 
transition. In order to handle the case properly, we need 
an intermediate bouncing phase B which smoothly con¬ 
nects the two phases / and II. We consider the collapsing 
(J) and expanding (II) phases smoothly connected by a 
nonsingular bouncing phase (B). Assuming the curva¬ 
ture is not important in phases I and II near the bounce, 
and assuming wi and wu for the two phases we have 

ai(t) = a I0 [~(t - f/)] 3 d+”/>, (47) 

an(t ) = a H0 (t - t n ) 3 < 1 +’"ju> . (48) 


The coefficients should be determined by matching a and 
a at the transitions + and t 2 - 
In the expanding phase II, the d -modes of 4>, tp x , tps, 
ip v and tp K in eqs. (43-45|) decay away whereas C-modes 
remain constant and have the roles of the relatively grow¬ 
ing modes. Whereas, as t —> ti in the contracting phase 
although the C-modes remain constant, the d-modes di¬ 
verge (for — 1 < wi < 1) as 


1 1-wi 

$ OC p v oops OC Pk OC OC \t-tl\ 1+W I 

H _5+3jvj_ 

<p x oc — OC \t — tj\ 3 + 3 “/ . 


(49) 


For wi = 0, the d-mode of $ vanishes exactly, and the 
d-mode of ip v vanishes in near flat situation. For wi = 1, 
the d-modes of $, <p v , ps and p K show In \ t — ti\ diver¬ 
gence, instead. A complete set of solutions in several 
different gauge conditions is presented in Tables 2-5 of 
p5[ ; although the solutions were derived in the expanding 
phase, the same solutions remain valid in the collapsing 
phase with the time replaced by its absolute value. 

A simple example of the bounce is the case with I\ > 0 
and a positive A |mJ 


as{t) = -\/3A7Acosh(y / A/3t). (50) 

Evolution of the scale factor is plotted in Fig. ||. Either 
for the vanishing fluid with pure A, or for a A-type fluid, 
we have p + p = 0. If we have p + p = 0 strictly the 
basic set of perturbation equations becomes trivial and 
we cannot determine the perturbations properly, i.e., we 
do not have meaningful perturbations. 


(5 + 3w/)(5 + 3 wu) Hi 


(46) 
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FIG. 2. Evolution of the scale factor (line). The dotted 
lines indicate the extensions of the collapsing, bouncing and 
expanding phases without matchings. The bouncing phase 
lasts for |t_e| < 5 with A = 1 and K = 1. As an example, we 
take wi = 0 and wn = |. 


In the following we consider the perturbation evolu¬ 
tions in three other examples of the bouncing phase. The 
first two models rely on the fluid/field which give effec¬ 
tively w < — | equation of state during the bounce. To 
have the bounce in such models the positive curvature 
should have significant role during the bounce, thus not 
suitable for the bouncing model assumed in §1V. In ad¬ 
dition, as the background curvature becomes important, 
all the perturbation scales go through the small-scale 
regime where the pressure gradient term becomes impor¬ 
tant. The third model relies on a presence of an exotic 
matter which gives a negative contribution to the total 
energy density. In this case we have a bounce without re¬ 
sorting to the positive curvature, thus the scales remain 
in the large-scale during the bounce and the model suits 
the requirement of the bounce assumed in §IV. 


A. Bounce with a w = — ~ fluid 


For an ideal fluid with w = constant, K > 0 and 
A = 0, from the Friedmann equation we have H = 0 at 
a(t*) = (2co/K) 1 ^ 1+3w ^ where Co = (47rG/3)/xa 3 ( 1+lu ^. 
We can show that a(t*) is a maximum for w > — 1, and 
is a minimum for w < — 4. As a simple example which 
gives a bounce we consider w = — | case, |35). Although 
it is uncertain whether it is appropriate to consider an 
ideal fluid for w < 0 case, we will take the ideal fluid 
assumption, see a cautionary remark in §VII of |J. We 
will find a fundamentally different result in a more real¬ 
istic (in the sense that we have the concrete action and 
equation of motion) case based on a scalar field, see § VB. 


Equation (||) gives an exact solution 

a = ( c o/2)(t 2 + K/cq). (51) 

We have c 2 = — | and 


X 

X 


cl 3f 4 + {K/c 2 ) 2 
2 a 2 t 2 




3 K/cl). (52) 


The pressure terms become important compared with the 
gravity near the bounce. Thus, in order to follow the 
evolutions we need to handle perturbations based on the 
full equations. Equation (|l4]) gives 


4t . 8 k 2 - 3 K 

^ + t 2 + K/cf x ~ 3 f 0 {t 2 + K/c 2 ) 2 * x ~ °- 

Ignoring the k 2 term the exact solutions in eqs. 
can be integrated, and for ip x we have 


(53) 


( 160 ) 


It 4 + 6 (K/cl)t 2 - 3(K/cq) 2 
3 ( t 2 + K/c 2 ) 2 


+d 


it/c 0 

(t 2 + K/ Cg) 2 ' 


(54) 


Since the k 2 terms become negligible away from the 
bounce we can use these solutions as the proper initial 
conditions for the C- and d-modes. A typical evolution 
is presented in Fig. [| 



FIG. 3. Evolutions of 4> for the G-mode (blue, line) 
and d-mode (red, long-dashed line) initial conditions, and tp x 
for the G-mode (cyan, dotted line) and d-mode (magenta, 
dot-short-dashed line) initial conditions. We take n = 10. 
The pressure terms become important in |t| < .4 for and 
in |t| < .2 for <E>. Notice that the G-mode of 4> changes sign 
twice, whereas the d-mode changes once. The sign changes in 
the expanding phase occur at the same time. 
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As we have we anticipate an exponential 

growth/decay of the perturbation while the pressure gra¬ 
dient term becomes important. In Fig. || both the C- and 
d-modes in the contracting phase become the (relatively 
growing) C-mode in the expanding phase. As the large- 
scale conditions are violated both the C- and d-modes 
will be dominated by the exponentially growing mode in 
the small-scale limit, and eventually we cannot trace the 
C- and d-modes in the expanding phase to the ones in 
the contracting phase. 


B. Bounce with a massive scalar field 


We consider a bouncing model based on a massive min¬ 
imally coupled scalar field with a positive curvature |3ti[ . 
Results up to eq. © in §|n] remain valid for the field 
with a prescription mentioned below eq. ( 0 ). The back¬ 
ground equations are presented in eq. 0 and below it 


with HQ = ^(/> 2 + m 2 (j) 2 ), etc. Equation (14) gives 


<Px' 


m 2 d>\ 

7 H + 

m 2 d> 2 „ rr m 2 d> 


k 2 - 8 K 


<Px = 0, (55) 


where M 2 t = 1/(87 tG). Once we have ip x , the rest of per¬ 
turbations 4>, tp v , (ps and if K follow from eqs. fl.|^. [To| . |TT| ) . 
We have 



FIG. 4. Evolutions of ip x (red, long-dashed line), <t> (blue, 
line), <p v (cyan, dotted line), and ip K (green, dot-short-dashed 
line). We omit i ps which has a larger amplitude and ap¬ 
proaches <f> in the large-scale. We take n = 10. t = 0 corre¬ 
sponds to the minimum of the bounce, and the pressure dom¬ 
inates over the gravity till t ~ 1. We can show that <t>, <p v , 
ip K and ips stay constant in the large-scale, whereas <p x still 
adjusting its value even in the large-scale as the background 
equation of state effectively changes from w<j> — — 1 while <f> 
is rolling to = 0 as <j> starts oscillating which occurs for 
t > 11. 


9 k 2 2m 2 d> K 

a 2 H<j) a 2 

y 2m 2 cj> I\ 2 

- =-:-9 + m 

y H(j> a 2 


k 2 


a 


2 ’ 


3 I< 

4^2 


(56) 

7cj) 2 + 25m 2 /) 2 
24 


+2m 


2^ 

$ 



4 H 


(57) 


We can show that near the bouncing era the pressure 
term in eq. (^b|) dominates over the gravity term in eq. 
([S’?]); this is true even for k 2 = 0. Thus, near the bounce 
the large-scale condition is violated, and with the pos¬ 
itive sign in front of k 2 term in eq. (^6|) we can show 
that perturbations show oscillatory behavior while in the 
small scale, see Fig. |1|. Although the first term in the 
RHS of eq. (56) diverges near the bounce, the same term 
appears in the gravity part in eq. ([57]) as well. Due to 
the positive sign in the second term we expect oscillatory 
instability as the term dominates the gravity part. Near 
the bounce we have H ~ 0, thus ~ constant, and 


a ~ J3M 2 l K/n (j> cosh J /^/(3 M 2 t )t 


(58) 


In this model, during the bounce all scales reach the 
small-scale where we cannot apply our large-scale solu¬ 
tions. In Fig. 0 we used an arbitrary initial condition at 
the minimum of the bounce (t = 0), and as the scale be¬ 
comes large-scale we have only the C-mode because the 
d-mode in expanding phase is decaying (thus transient) 
and yields to the relatively growing C-mode within a few 
expansion time scale. Thus, as in the previous example 
based on the w = — | fluid the C- and d-modes during the 
collapsing phase are mixed up while in the small-scale, 
and we cannot trace the C- and d-modes in the expanding 
phase to the ones in the collapsing phase. One important 
difference of the scalar field compared with the w = — | 
fluid is that, while fluctuation of the fluid shows expo¬ 
nential instability due to the negative c? term, the field 
fluctuation shows oscillatory behavior |3j| . Such a differ¬ 
ence comes from the presence of a nonvanishing entropic 
perturbation term e in the case of a scalar field, see below 
eq. (|). 

In later expanding phase as <fi starts oscillating near 
potential minimum the background model enters an era 
with effectively = 0 (dust) equation of state. As 
/> starts oscillating we cannot solve eq. (|5^) directly. 
Instead, we can analytically handle the situation using 
proper time averaging over the coherent oscillations of 
the background and the perturbed field |Q. In [|0| it 


















was shown that while the background enters the dust 
era, the perturbations also behave like a cold dark matter 
even in the large-scale limit. 


C. Bounce model with fi = fi m — fix 


The positively curved FLRW world model with only 
the radiation and matter does not allow bouncing after 
the big crunch. If the physical state near the big crunch 
allows a presence of additional matter X with its effective 
energy density behaving as —px(t) = — pxoo-~ 3( ' 1+wx ' 
and wx > 5 , we could have a smooth and non-singular 
bouncing phase; this is a generalized case of a ‘desparate’ 
example mentioned in p368 of j||. Thus, for the bounce 
only, we do not even need the positive curvature in the 
background world model. 

As a toy model allowing such a smooth and non¬ 
singular transition with the relevant scale satisfying the 
large-scale condition we consider a case with the pressure¬ 
less matter and the exotic matter with u>x = f • Thus, 
we consider a model with the dust and the radition with 
a negative sign in the radiation component. 

„ 2 87tG . K 

H = —X-(Vm - Px) - 2 > ( 59 ) 

where fi m oc a -3 and fix oc a -4 . Certainly, this is not 
a realistic model for the bounce because we need to as¬ 
sume that there is no conventional radiation component 
present; to be more realistic, the sum of the radiation 
and the X-matter should give a net negative contribu¬ 
tion. Later we will show, however, that this toy model 
captures the basic physics of more realistic situations. 

For a positive curvature, K > 0, eq. (^) gives an 
exact solution 


^ - T r Cn 

K 


1 - \J\- 2(c x /cl l )K cos {VKr]) 


(60) 


where cx = ( 47 rG/ 3 )/xxa 4 ; Cx/c ^ is dimensionless. We 
have normalized the time axis so that a has minimum 
values at rf = 2fnr/\/H with n an integer number. For 
vanishing X-component, cx = 0, we recover the solution 
in eq. ([32|). With the X-component the model shows a 
cyclic behavior. The basic picture of the cyclic bounces 
remains valid in more realistic situations with wx > y • 
The K term becomes important near a max , and near 
< 2 min we have fix — Pm- The curvature term is negligible 
near the bounce, thus allowing existence of the large-scale 
where we could ignore the Laplacian term coming from 
the pressure gradient. If we do not need a cyclic behavior 
we can actually take a flat model with K = 0. Assuming 
negligible K term we have 


+ 2c x A4), 


(61) 


which is an exact solution of eq. (|59|) for K = 0. 


It is not entirely clear how to handle the perturbation 
of the exotic component X which is introduced to have a 
bounce in the background without resorting to the posi¬ 
tive curvature. If such an exotic state of matter can be 
modelled by using a held or modifying terms in the grav¬ 
ity theory (coming from quantum corrections or higher 
dimensions, for example), we need the correct forms of 
the equation of motion or the modified action to handle 
the perturbations properly. In the present situation lack¬ 
ing the concrete model for the exotic matter, we will take 
a phenomenological approach based on a fluid approxi¬ 
mation. 

Let us assume, except for its negative contribution to 
the density, the fluid X behaves as an ordinary ideal fluid 
with an equation of state wx- In such a case, near the 
bounce we need to consider two-component system with 
m and X. [The other possibility is to derive solutions 
in the to- and X-dominated eras separately and to con¬ 
nect by using the matching conditions. As the X-fluid 
cannot dominate the total fluid (we have fi m > fix with 
the equality holding at the bounce) such a situation is 
forbidden.] 

We use the conventional decomposition of the sys¬ 
tem into the adiabatic (curvature) perturbation, char¬ 
acterized by ip x (or <f>), and the relative (often called 
the isocurvature) perturbation defined as S = S 12 = 
5i/(1 + wi) — 52/(1- + W 2 ), The basic equations in the 
two-component system become, see eqs. (23,35,46,57) in 

PI 


P+P 

H 


H 2 


i{fi+p) 


(>*) 


+ c l^x = -4tt Ge 


= -47rG 


{pi +P1X/P2 + P2), 2 2 


P + P 
-,2 ' 


( c (l) c (2))*^) (62) 


S + H (2-3c 2 ) S + cl-^S 

k 2 (k 2 — 3 Ii) c 2 1 } -c 2 2) 
a 4 47 rG(/i + p) 




(63) 


where c 2 ^ =pi/fa, and 


2 _ C( 2) (pi +Pl) +C 2 1 } (/X 2 +p 2 ) 
p + p 


(64) 


In our case 1 = to and 2 = X, thus, wi = c 2 ^ = 0, 
W 2 = cf 0) = p + p = Pm — | px, etc. In the large- 
scale limit the curvature mode can be sourced by the 
isocurvature one, see eq. (|62|), whereas the isocurvature 
mode decouples from the curvature mode in general [ [Ill 1 
see eq. ([63]). [The situation is different in the case of 
multiple number of scalar held: in such a case, effectively, 
the RHS of eq. (^3|) has k 2 /a 2 factor instead of fc 4 /a 4 , 
thus the isocurvature modes are less decoupled from the 
adiabatic one, see US-] 

Let us assume an adiabatic initial condition, thus set¬ 
ting S = 0 at early era in the contracting phase, for 
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simplicity. More precisely, we are assuming S -C ip x in 
the initial epoch which is natural because it means we are 
assuming no significant fluctuations in the X-component 
in the early matter dominated era. Since the curvature 
mode does not source the isocurvature mode in the large- 
scale, the isocurvature mode will remain small. In such 
a case the RHS of eq. vanishes, and the curvature 
equations in the single component situation, eq. (14), 
remain valid without any change. Thus, for scales satis¬ 
fying the large-scale limit we have the same solutions in 
eqs. (^6 17) remaining valid. The solutions become: 


VI. DISCUSSIONS 




1 Px = 


1 + + 

2 c A' 


V 

5ci 


124 


4c 2 y 


4 ca" z -'-a \ ^ j 

wprv - — ) rj d V 


Q r 6 

40c 3 x 


V 6 c 


c 6 - 
c x 


1+ 2cA ?r 


(65) 


For <p x the contribution from the lower bound of integra¬ 
tion of the C-mode is absorbed to the d-mode. In the 
matter-dominated eras, |??| 2c x /c'^ n , we have 


135 clL 


ip x = -C + 5 d. 

0 


( 66 ) 


Near the big bang/crunch, the solution for ip x coincides 
with the one for a pressureless medium considered in eq. 
(p8|). For cx/c^n —> 0, $ also coincides with the one 
known in the pressureless medium. Near the bounce, 
\rj\ <C 2ca/c^, we have 

$ = C + k 2 ^frfd, tp x = C + ?fnd, (67) 
9 c x c x 

which are regular and finite. 

Since the present bounce model allows the scales to 
stay in the large-scale limit during the transition, it can 
be considered as a concrete model of the smooth and non¬ 
singular bounce assumed in § IV. Indeed, the curvature 


term is negligible near the bounce as it was the case near 
the big crunch/bang injfy]. In the cx/c^ —> 0 limit, eq. 
(|65|) reduces to eq. (|66|) which also coincides with the 


known solution considered in §IV. 

Apparently, we can also make a more realistic model 
with the radiation, matter and X where Wx > -- ■ 
We note that eqs. (|62|,|63|) remain valid for any two- 
component system of matter perturbations. Even in such 
a case the X fluid can cause a smooth and non-singular 
bounce and the curvature term has negligible role near 
the bounce. Thus, essentially the same conclusions (e.g., 
the C-mode feeding the C-mode) remain valid. We have 
considered a simple toy model only because it allows an¬ 
alytic handling of the background and the perturbations, 
thus showing the situation explicitly. 


In § IV we have shown that the perturbation in a posi¬ 


tively curved FLRW model filled with a pressureless mat¬ 
ter is described by the conservation of <F. Assuming a 
transition of big crunch followed by a big bang in such 
a model, by using the known matching conditions we 
have shown that $ maintains the same value even after 
the transition. Using the matching conditions we also 
have shown that the diverging solution in the contracting 
phase is matched to the decaying solution in the subse¬ 
quent expanding phase, whereas the other solution which 
stays constant during the contracting phase is matched 
into the same constant solution in the expanding phase. 
That constant mode is characterized by 4); the other solu¬ 
tion of $ which decays in expanding phase is higher-order 
in k 2 compared with the one of tp x and vanishes for the 
vanishing background pressure. 

In order to confirm these results based on the match¬ 
ing at singular bounce, in §0 we have considered three 
different non-singular and smooth bounc e mo dels. For 
the bounce models based on the fluid (§VA) and the 


massive scalar field (§VB) in the positively curved back¬ 
ground, the role of background curvature is important 
to make the bounce. In such cases, all the perturba¬ 
tion scales come inside the sound-horizon near bounce, 
and the large-scale conditions are violated. As the pres¬ 
sure gradient terms become important, perturbations in 
the fluid model show exponential instability, whereas the 
ones in the massive field model show oscillatory behav¬ 
iors. For both situations the two independent pertur¬ 
bation modes in the large-scale limit during the collaps¬ 
ing phase got mixed up with the two independent modes 
in the small-scale during the bouncing phase. Thus, we 
cannot trace the two independent solutions (C- and d- 
modes) in the expanding phase to the ones in the con- 
tracti ng p hase. 

In § |V C] we considered a bouncing model based on an 
exotic matter with negative contribution to the total en¬ 
ergy density. In such a case the positive curvature is 
not important during the bounce. Thus, we could have 
the relevant scales remaining in the large-scale limit, and 
could apply the general large-scale solutions. In this case, 
however, we have to handle the perturbation of the ex¬ 
otic matter in addition to the ordinary one simultane¬ 
ously. By considering the adiabatic initial condition we 
have shown that the same curvature perturbation equa¬ 
tion known in the single-component situation remains 
valid, thus the known large-scale solutions are valid as 
well throughout the bounce. Therefore, this third-type 
of bounce model can be regarded as an example of the 
smooth and non-singular bouncing assumed in § IV. As 
an analytically manageable concrete example, we have 
considered a simple toy model with a dust and an exotic 
matter with the radiation-like equation of state. Even in 
more general situations of the equation of states before 
and after the bounce, similar analyses can be made which 
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show that the C-mode in the expanding phase is affected 
only by the C-mode in the contracting phase, thus the 
growing mode in the contracting phase decays away as 
the world model enters expanding phase. 

Our analyses are based on two assumptions: (i) the 
contracting phase is converted into the expanding one 
by a smooth and non-singular bounce, and (ii) the lin¬ 
ear perturbation theory holds during the evolution. The 
large-scale evolution can be characterized by the conser¬ 
vation of <f>. We have shown that the C-mode of $, which 
is the proper growing mode in expanding phase, is sim¬ 
ply conserved during the evolution and through bounces. 
The results are true as long as the two assumptions made 
above are valid, and in addition, if the large-scale condi¬ 


tion is met during the transition as considered in §V_C 


In §111 we showed that the three dimensionless mea¬ 
sures, the intrinsic curvatures (p and C^g), the tracefree 


part of the extrinsic curvature (a / H ), and the Weyl cur¬ 
vature ( E/R ), diverge at singularity for — 1 < w < 1. 
Thus, for — 1 < w < 1, the spacetime perturbations be¬ 
come singular as the background approaches the singu¬ 
larity. An ambiguity remains for w = 1 case because 
although p and C^ diverge logarithmically, <j/H and 
E/R remain finite. These apply to all perturbation types 
and for all gauge conditions we have considered. Behav¬ 
iors of the other variables (the perturbed lapse function 
a, a dimensionless measure of the perturbed expansion 
k/H, relative density perturbation S , etc.) depend more 
strongly on the gauge conditions, see Tables 2 and 4 of 
p5[ . Thus, these variables apparently have less physi¬ 
cal significance to characterize the spacetime fluctuations 
compared with the other three measures whose behav¬ 
iors are gauge independent at least in the pool of gauge 
conditions we have investigated. Do above results imply 
diverging spacetime fluctuations for — 1 < w < 1, and 
regular ones for w = 1? In Table 4 of (^5| we find that in 
no gauge condition we have all the perturbations remain 
finite for — 1 < w < 1. 

The authors of [Q argued that as the model goes 
through a singular bounce the perturbation becomes non¬ 
linear. We have shown that if the fluctuations survive the 
bounce as linear ones, the diverging mode in the contract¬ 
ing phase should be matched to the decaying one in the 
expanding phase. Lyth in ||44|| made the following sim¬ 
ple and powerful argument. As we have under the gauge 
transformation ( x a = x a + with £f = af v ) 

p = p — Hf t , S = 6 + 3H(1 + , (68) 


if p diverges while 6 remains finite, or vice versa, in any 
single gauge condition [this is the case for — 1 < w < 
1, see eqs. (g|49) for p v and <5„] no temporal gauge- 
transformation can be found which makes both p and 
S finite. Therefore, for — 1 < w < 1 we encounter the d- 
mode perturbations of Friedmann world model becoming 
singular near big-crunch in one form or another in all 
gauge conditions. 


We note that <f>, which becomes p v for K = 0, sim¬ 
ply stays constant in a pressureless medium, thus its 
magnitude cannot characterize the breakdown of linear¬ 
ity of the perturbation. As we have from eqs. (f^|. [i7|Jd^ ) 
p v = C + (1/3)(1 — cosrj)p x where we set K = 1, p v 
itself could diverge near singularity. From eq. (|l^), near 
big crunch in the pressureless medium we have p K ~ p$ 
where ps, given in eq. (39), has diverging part. Thus, 
near big crunch the diverging modes behave as 


Pv oc p s oc p K oc \rj\ / p x oc \rj\ 


1-5 


(69) 


whereas <f> has no diverging mode in the pressureless case. 
For the situation with general w, see eq. (|4^). Bardeen 
has argued that the behavior of p x “overstates the phys¬ 
ical strength of the singularity”, see below eq. (5.12) in 


At the singular big crunch, we certainly have d-modes 
of many perturbation variables unambiguously becom¬ 
ing singular for —1 < w < 1, see Tables 2-4 in 125 . Do 


large amplitudes of perturbations imply breakdown of 
linear theory? Due to the gauge dependence of relativis¬ 
tic perturbation, large (larger than unity, say) amplitudes 
of some dimensionless measures of gauge-invariant per¬ 
turbation variables do not guarantee the breakdown of 
the linear theory. However, what Lyth Q] has shown is 
that in the collapsing phase we could encounter situations 
where the amplitudes of perturbation variables become 
large in one form or the other in all gauge conditions. 
Lyth has argued this as the violation of the necessary 
condition of the linear perturbation theory. 

In our models avoiding the singularity by a smooth 
and nonsingular bounce it is likely that certain scales can 
safely go through the bounce retaining their linear na¬ 
ture. As we have assumed the linearity of perturbations, 
our analyses and results are applicable to such scales only. 
In §[V and V C we have shown that the diverging solu¬ 
tions in the contracting phase in eq. (|6^) affect only the 
decaying, thus transient, solutions in the subsequent ex¬ 
panding phase. In such a scenario, however, one could 
anticipate large (compared with the C-mode) amount of 
the decaying ( d ) mode present in the early big-bang phase 
for a while which is the remnant from the preceding phase 
before the big-bang. 

In a recently proposed ekpyrotic scenario based on col¬ 
liding branes Q it was argued that the final scalar-type 
perturbation is scale invariant |4(|. In [[idj it was shown 
that the generated scale-invariant spectrum in the zero- 
shear gauge during the collapsing phase should be iden¬ 
tified as the d-mode, thus after the bounce we have a 
different power-spectrum |48[ ]. Our results in this paper 
confirm that while the large-scale condition is met during 
the (smooth and non-singular) transition the d-mode in 
the contracting phase does not affect the (properly grow¬ 
ing) C-mode in the expanding phase. The background 
curvature is flat in the ekpyrotic scenario and the scale 
remains in the large-scale during the bounce. However, 
since the bounce of the ekpyrotic scenario goes through 
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a singularity the author of jlij has argued that one can¬ 
not rely on the linear analyses as the model approaches 
the singularity. Thus, either the final spectrum is not 
scale-invariant (which is the case if the linear perturba¬ 
tion survives) or the issue should be handled in the future 
string theory context with a concrete mechanism for the 
bounce. 
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